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■ Abstract 

, ^ . For any graph consisting of k vertices and m edges we construct an ensemble of random pure 

' quantum states which describe a system composed of 2m subsystems. Each edge of the graph 

represents a bi-partite, maximally entangled state. Each vertex represents a random unitary 
matrix generated according to the Haar measure, which describes the coupling between subsystems. 
' Dividing all subsystems into two parts, one may study entanglement with respect to this partition. 

A general technique to derive an expression for the average entanglement entropy of random pure 
states associated to a given graph is presented. Our technique relies on Weingarten calculus and 
flow problems. We analyze statistical properties of spectra of such random density matrices and 
^ ^1 show for which cases they are described by the free Poissonian (Marchenko-Pastur) distribution. 

We derive a discrete family of generalized, Fuss-Catalan distributions and explicitly construct 
, graphs which lead to ensembles of random states characterized by these novel distributions of 

■ eigenvalues. 
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1 Introduction 



The phenomenon of quantum entanglement in physical systems remains a subject of considerable 
X/^ ' scientific interest. The case of entanglement in bi-partite systems is relatively well understood [T], 

. but the case of systems consisting of several subsystems is much more demanding and complicated. 

' For instance, the measures characterizing quantitatively quantum entanglement worked out for the 

simplest case of two subsystems are often not capable to describe the complexity of multi-partite 
CO ' entanglement. On the other hand, the measures applicable in this case, like the geometric measure 

. of entanglement [4j related to the minimal distance of the analyzed pure state to the set of separable 

states, are in general not easy to compute. 

Not being in position to characterize entanglement of a concrete state of a composed quantum 
system, one often tries to describe properties of a 'generic' quantum state. To this end one constructs 
ensembles of quantum states and computes interesting quantities averaged over the entire ensemble. 
' In the simplest case of a system consisting of two subsystems one considers an ensemble of ran- 

dom pure states distributed uniformly with respect to the natural, unitarily invariant (Fubini- Study) 
measure. Quantum entanglement with respect to the partition into these subsystems called A and B 
can be described by its entanglement entropy, E{(p) = S{TrB \(p){(p\). Here S{a) = — log Tr u log cr 
denotes the von Neumann entropy of the state a obtained by the partial trace of the state \lp){lp\. As 
the initial state is pure, the resulting entropy does not depend on the subsystem, with respect to which 
the partial trace is taken. The entropy is bounded from above by the log of the smaller dimensionality, 
5'max — log dm, whcrc dm = mm{dA, ds} ■ If this bound is saturated, the initial pure state is called 
maximally entangled. 

The mean entanglement entropy {E{ip))^, averaged over the ensemble of random pure states de- 
pends only on the dimensionality of both subsystems. An explicit formula for the mean entropy of a 
subsystem was first conjectured by Page [SI and later proved in |^. Due to the measure concentration 
phenomenon the entanglement entropy of a generic pure state is close to the mean value (E) . Concrete 
bounds relying on the Levy's lemma [7] specify the probability to find a random state with entropy of 
entanglement smaller than (E) by e. 



1 



Note that the Fubini-Study measure on the space of pure states of the entire system followed by the 
partial trace over subsystem B induces a probability measure in the space of mixed quantum states 
describing the subsystem A. These induced measures [8l[9] are parameterized by the dimensionality of 
the auxiliary subsystem B. 

Analyzing an ensemble of the reduced random states cr, Page derived the probability distribution 
for its eigenvalues. Up to an overall normalization constant this problem is equivalent to finding the 
density of the spectrum of the Wishart matrices W = GG* (where G is a rectangular complex random 
matrix of the Ginibre ensemble) studied earlier by Marchenko and Pastur [10^- This very distribution 
depending only on the ratio c — ds /(Ia of dimensions of both subsystems, is also called free Poissonian, 
as it corresponds to the free convolution of random matrices [11(112) . 

In the general case of pure random quantum states describing multipartite systems one usually 
studies entanglement with respect to various cuts of the entire system into two parts |13ffT5] . so the 
standard measures of the bi-partite entanglement can be applied. Statistical properties of random pure 
states change if one breaks the overall unitary invariance, characteristic to the Fubini-Study measure, 
and introduces some structure into the multi-partite system. In particular, for certain numbers of 
subsystems there exist perfect maximally multipartite entangled states (MMES), such that they are 
maximally entangled with respect to any bi-partition [14,16 . 

Random quantum states are useful to tackle various problems of theoretical physics. A key con- 
jecture of the theory of quantum information processing concerning the additivity of minimal output 
entropy was recently shown to be false |17II18| . It is worth to emphasize that the original reasoning was 
not constructive but it was based on relations between the average quantities computed with respect 
to an ensemble of quantum states. Entanglement between random states is also interesting in systems 
motivated by condensed matter physics 19,20]. For a discussion of a relationship between quantum 
criticality, Anderson transition and the average entanglement entropy see the recent review |21| and 
references therein. 

The theory of random pure states proved to be useful to analyze the information flow and the 
entropy of black holes. The model of Page [22] consisted of two subsystems, the dimensions of which 
served as parameters of the model. In later models one studied pure states of a system composed 
of four parties 23 , 24 . The system contains two pairs of maximally entangled states which relate 
subsystems A, A' and B, B' respectively. The subsystems A and B are directly coupled together, and 
the action of an unknown Hamiltonian is mimicked by a random unitary matrix distributed according 
to the Haar measure. Such a system can symbolically be depicted by a graph containing two edges, 
which represent two maximally entangled states and a vertex denoting a random unitary matrix. This 
very situation is shown in Fig. [5]d, where a slightly different notation is used. 

The main aim of this paper is to extend this construction for an arbitrary (undirected) graph and 
to investigate properties of the resulting ensembles of multi-partite random quantum pure states. Any 
graph with m edges will be associated with a quantum system consisting of 2m subsystems. Each edge 
corresponds to a maximally entangled state while each vertex represents a random unitary matrix. 
The dimensions of two subsystems linked by a given edge of the graph are assumed to be equal, but 
besides this constraint the dimensions of all the subspaces can be treated as free parameters of the 
model. As the physical interaction is modeled by a set of random unitary matrices each graph defines 
an entire ensemble, we call them random graphs states. 

Note that such an ensemble of quantum random states associated with a given graph differs form the 
deterministic construction of graph states introduced by Hein et al. [25] . Furthermore, our approach is 
not related to quantum graphs reviewed by Gnutzmann and Smilansky [37], which describe quantum 
particles (or waves) traveling along a graph. A graphic representation of ensembles of random states 
analyzed in this work looks slightly similar to the one used to define the Projected Entangled Pair 
States (PEPS) formalism, capable to describe complex many-body systems [28|I29] . While in the latter 
setup the subsystems entangled with auxiliary systems are coupled together by a projection on a low 
dimensional subspace, in our approach such a coupling is described by a random unitary matrix. 

Any graph defines the topology of couplings between the physical subsystems. Choosing a certain 
set of the subsystems one can average the pure quantum state {ip) over the remaining subsystems. Tech- 
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nically one performs the partial trace the auxiliary subsystems, which leads to a reduced state which 
generically is not pure. Hence selecting a concrete graph and specifying the subsystems pertaining to 
the environment we define an ensemble of random mixed states. 

In this work we develop general techniques suitable to describe spectral properties of random 
density matrices associated with a graph. Our technique relies on Weingarten calculus [30, 31 , and 
the observation that its asymptotics boil down to a min-flow problem [32 . In the asymptotic limit, as 
the dimension of the quantum states goes to infinity, statistical properties the spectrum of a random 
state can be established, so explicit formulae for its purity and entropy are derived. 

We identify these ensembles of random states for which spectra are described by the free Poisson 
(Marchenko-Pastur) distribution. In some other cases the spectra are described by a discrete family 
of probability distributions 7r(*\ which can be be considered as generalizations of the free Poissonian 
distribution. We call them Fuss- Catalan distributions (FC), since their moments are related to the 
Fuss-Catalan numbers, known in the combinatorics and free probability calculus [331134] . 

The Fuss-Catalan probability distributions are also closely related with properties of products of 
non-hcrmitian random matrices. In general, products of random matrices are often studied in context 
of various problems of statistical physics |35j . Recent studies on products of Ginibre matrices concern 
multiplicative diffusion processes [36], rectangular correlation matrices used in macroeconomic time 
series [37j and lattice gauge field theories [3B]. Spectral properties of a product of random Ginibre 
matrices were analyzed in a recent paper of by Burda et al. [39]. 

In the simplest case s — 2 random matrices described by the distribution tt'^^' have the structure 
G2GiGlG'2 where Gi and G2 are independent random Ginibre matrices. In the general case of an 
arbitrary integer s the random states are proportional to the Wishart matrix W = GG* , constructed 
out of a product of s independent Ginibre matrices, G — Y[i=i ^i- 

To analyze statistical properties of random states whose definition is based on random unitary 
matrices, one needs to perform averages over the group of unitary matrices. Instead of using explicit 
formulae derived by Mello [40] we found it more convenient to use the Weingarten calculus (see e.g. 
[30ll41j). In order to evaluate expectation values for random tensors we are going to use a diagramatic 
approach and a graphical calculus recently developed in |31] . 

This paper is organized as follows: ensembles of pure states corresponding to classical graphs are 
defined in section [2| In section [3l we discuss mixed states which arise by taking the partial trace 
over a specified set of subsystems. In section |4] graphical and combinatorial tools used to perform 
integration over the unitary group are described. The main result of this work — a general technique 
to compute the moments of the spectrum of ensembles of mixed states obtained by partial trace — 
is presented in section [S] In section [5] we analyze ensembles associated to some particular graphs 
(star graphs, cycle graphs) and find ensembles of mixed states characterized by Marchenko-Pastur and 
Fuss-Catalan distributions. Some other graphs which lead to other, "exotic" distributions are also 
provided. The paper is concluded in section 7 in which we summarize results obtained and present a 
list of open questions. 

2 Classical graphs, quantum pure states and random matrices 

In this section we describe a family of quantum pure states that are associated in a natural way with 
undirected graphs. The underlying idea is that vertices of graphs correspond to quantum systems and 
edges describe the entanglement between the systems. A particular feature of the states we consider 
is that the entanglement described by the edges will be maximal. 

Consider an undirected graph T consisting of m edges (or bonds) i3i , . . . , _B,„ and k vertices 
Vi,...Vk- We allow multiple edges between two vertices, as well as vertex loops. Note that the 
graphs considered here are not metric, so we do not discuss the lengths of the edges. Let bi denote 
the degree of the vertex Vi defined as the number of edges attached to this vertex (each loop counts 
twice). As an example, for the graph F in Figure [U we have 61 = 1, 62 = 2, 63 = 3. 

We are going to analyze quantum states belonging to the Hilbert space with the n — 2m-fold 
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tensor product structure, % = Hi 'Si ■ ■ ■ ^ 'H2m- The dimension of the subspace Hi is diN , where the 
parameters di > Q are fixed while the number N is arbitrary and we will eventually discuss the limit 
— > oo. The total dimension of space Ti. reads thus D = dini'H = (^nj=i 

A graph T with m edges induces a tensor product structure of the 2m- fold tensor product space H. 
The vertices Vi of F induce a partition Ilvortox of the set [n] = {1, 2, . . . , 2m}, with blocks Ci, . . . , Cfc. 
The reader can think about the Hilbert spaces {T-lj \ j E Ci} living inside the vertex V/. For each 
vertex Vi, we introduce the Hilbert space Wi — Sj^CiT~l-j which describes the 5/ subsystems of the 
vertex Vi . 

Each edge of the graph represents the maximally entangled states between two Hilbert spaces Hi 
and Hj. More precisely, we shall ask that dim Hi = dimHj — diN and we put 

diN 

H) = ^T.\^.)®\h). (1) 

^ ' x=l 

where \ex) and \fx) are orthonormal bases for Hi and Hj. The actual choice of the two bases is not 
important in what follows, but we shall assume that these bases are fixed. 

In order to introduce an ensemble of random graph states belonging to H we are going to use k 
independent, random matrices Ui,U2-, ■ ■ ■ ,Uk, distributed according to the Haar measure. The matrix 
Ui corresponds to the l-th vertex of the graph, VJ. It acts on the space Wi (and thus on bi subspaces) 
and describes the coupling at a given vertex. To specify the subspaces coupled by a unitary matrix it 
will be also convenient to use an alternative notation and write Ui = Ui-^^,,,_i^ where p = bi denotes the 
degree of the vertex. See Figure [1] for a concrete example of a graph and a graph ensemble. 




(a) (b) 



Figure 1: (a) Graph consisting of 3 vertices and 3 bonds, one of which is connected to the same vertex 
so it forms a loop; (b) the corresponding ensemble of random pure states defined in a Hilbert space 
composed of 6 subspaces represented by dark dots. Each bond corresponds to a maximally entangled 
state spanned on two subspaces connected by the bond, while each vertex corresponds to a random 
unitary matrix. 

The simplest graph ensembles described by our model have only one edge. We consider two cases, 
one where the unique edge is a loop of one vertex, and the second one where the edge connects two 
vertices (see Figure [5]). 

In the first situation, the graph has a single vertex of degree two. This graph corresponds to the 
system consisting of 2m — 2 subsystems. The edge represents the maximally entangled state, but due 
to the global unitary matrix Ui ~ C/1,2, l^*) — C/i,2|*i'i2) ^ random state in Hi (8) H2- Since Ui^2 is 
generated according to the Haar measure, the random state j^") is distributed uniformly with respect 
to the Fubini-Study measure, as in the model of Page [5] , so in this case the bond forming the loop is 
in a sense redundant. 

In the linear case, we have again only one edge, m = 1, two vertices, fc = 2, and each vertex 
contains only one node. The graph is represented in Figure [2] (c) and (d), both in the simplified and in 
the standard form. The random quantum pure state associated to this graph is a generic maximally 
entangled state, 

dN ^ dN 

I*) = [Ul SU2]-^Y1 1^^) ^ 1/^) = E If^ie.) ^ \U2f^), (2) 
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where we set di = d2 = d and where Ui,U2 represent independent random unitary matrices. Here 
{ciji , {fi}i denote some fixed bases of C*^. 

8 T^.Qh, o o '^^^^ 

(a) (b) (c) (d) 

Figure 2: Graphs with one edge: a loop on one vertex, in simphfied notation (a) and in the standard 
notation (b), and two vertices connected by one edge, in simphfied notation (c) and in the standard 
notation (d). 

As a third and final example, the next simplest graph one could imagine has a linear shape, m = 2 
edges, k = 3 vertices Vi, V2, V3 and n — 2m — 4 nodes, as in Fig. [H 



Ti.1 H2 Ti-s H4 

(a) (b) 

Figure 3: A linear 2-edge graph, in the simplified notation (a) and in the standard notation (b). 

In general, any graph F is associated with a collection oi n — 2m subspaces "Hi, . . . ,"^2™, each 
endowed with some fixed orthonormal basis. We consider two partitions of the set [n] = {1, 2, . . . , n}. 
A pair partition Hedge with the property that V{i,j} G Hodge, di — dj, represents all bonds in the 
graph and thus the m entangled states. The second partition Hvertex consists of k blocks of size bi, 
which represent the vertices of the graph, and encode the random unitary operators Ui. For instance, 
in the case of n = 6 subspaces presented in Fig[T]the two partitions read Hedge = {{1, 2}, {3, 4}, {5, 6}} 
and Hvertex = {{!}, {2, 3}, {4, 5, 6}}. 

These two partitions allows us to introduce a general definition of an ensemble of graph random 
states 



OS) 

.cenvcrtox 



, {i,j}encdEc 




(3) 



where Uc G l4{Wc) — U{(^i£cT~l-i) are Haar independent random unitaries and ^^'^ the maximally 

entangled states H]) defined with respect to the fixed orthonormal bases of Hi and Hj. 



3 Marginals of random graph states 

To study non-local properties of the random graph state \'^) associated to a graph F it is useful to 
specify a partition of the set of all 2m subsystems into two groups, Htrace — {S, T}. Then it is possible 
to treat the multi-partite state j'l') as if it where bipartite and to analyze its entanglement with respect 
to this concrete partition. The entanglement of a bipartite pure state can be measured by the von 
Neumann entropy of the reduced state obtained from the initial projector by the partial trace. 

We will analyze the partial trace of the pure state over a subspace Ht defined by the subset 

T of the set [n]. Then S denotes the complementary subset, so the sum of their elements |r| + l^l is 
equal to n and the entire Hilbert space can be written as a tensor product, H = Ht Hs- Such a 
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partial trace will be briefly denoted by Try and the resulting density matrix is supported then by the 
subspace Hs, 



ps = TrTmm. (4) 

Graphically, partial traces are denoted at the graph by "crossing" the spaces T-Li which are being 
traced out. For example, in Figure SI the partition defining the bi-partite structure is Iltracc = {5* = 
{1,3, 5}, T = {2,4,6}}. In the diagram, the dots representing the subspace H2, ^4 and "He are 
"crossed", and the partial trace is taken on the tensor product of those spaces Tir = 'H2®'Ha®'Hq. 

Before exploring the quantitative entanglement of we can make some qualitative remarks 

right away. Looking at Figure |4] we see that the systems 'Hi and are not entangled before the 
random unitary transformation V4 = C/4,5.6 was applied. Indeed, before applying the random unitary 
matrices, the state of the system reads 

|^')-|<i>+)®|$+)®|$+). (5) 

The product structure of the state \^) is given by the supremum of the vertex and edge partition, 
Ilvcrtcx V Ilodgo. However, after the unitary transformations, subsystems "Hi and "Hs become entangled 
(generically) ; one might say that entanglement was broadcasted by the unitary matrix V4. 
In general, we define the partition Ilontangic by 

^ J « J and [i,i]n,.ertcx ^ (6) 

where [*,j]nvortox denotes the block of i and j in Ilvortcx. In other words, Ilentangic is obtained from 
Ilvertex in the following way: if a block of Ilvcrtox is entirely traced out, then replace it by singletons, 
otherwise keep it as it is. The block structure of ps is then given by the restriction of the partition 
Ilontangic V Hodge to the sct S. For the graph in FigurelH we have n — Q, Hodgo = {{1, 4}, {2, 5}, {3, 6}}, 
Hvortox = {{I}) {2}, {3}, {4, 5, 6}} and S = {1,3,5}. Since the unitary U^^^^ is not completely traced 
out, it is capable to broadcast the entanglement between subsystems 1 and 3. 



Hi 
^1® 



4( * 



^2 



^3® 



7^4 



^5 



Figure 4: Exemplary graph representing random states supported on n = 6 subspaces and defined by 
entangled states |<I'^4), |'I'25) ^^-^ |3>j^g) and a non-local random unitary matrix V4 ~ ?/4.5,6 of dimension 
d^d^d^N^ and local unitary matrices ¥±,¥2 and V3. Partial trace of this state over the subspace Ht 
defined by the set T — {2,4,6}, represented by crosses, provides the reduced state ps supported on 
subspaces corresponding to the set S = {1, 3, 5} and represented by full (non crossed) dots. 

We finish this section by a remark which can simplify the structure of graph state marginals in some 
cases. If a marginal of the pure state ([3]) is considered, the unitaries Uc corresponding to singletons 
C = {i} G Hvortox in equation ([3]) can be replaced by the identity operator if the pair j of i in Hodge is 
not a singleton in Hvertex- This follows from the invariance of the Haar measure on the unitary group 
and from the fact that one can extract an independent unitary from the block of j in Hvortox- 



4 Graphical and combinatorial tools for unitary integration 
4.1 Permutations, non-crossing partitions 

In this section, we introduce notation that will be used in the entire paper and which may be non- 
standard. For an integer n, let [n] denote a set of n elements {1, 2, . . . , n}. For a collection of numbers 
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{di}i, Hilbert spaces {Hi}i and for a set of indices /, we are going to use the notation Hi — ^ieiHi 
and di = Oie/'^*- 

We denote by Sq the group of permutations on q elements. For a permutation a, we call \a\ its 
length, that is the minimum number of transpositions necessary to obtain a and #(t the number of 
disjoint cycles of a. One has \a\ + #(t = q for all permutations a £ Sq. Note that the absolute value 
is also used in this work to denote the cardinality of a set \A\. 

We consider the lattice Vq of partitions of [g] = {1, . . . ,q}, and say that 11 ^ 11' if, for any block V 
of n there exists a block V' of H' that contains V, V CV' . For any partitions 11 and 11' define 11 V 11' 
(resp. n A n') to be the least upper bound (resp. the greatest lower bound) of 11 and H' with respect 
to the previously defined partial order, and also let Ig = {{1, . . . , q}} (resp. 0^ = {{1}, . . . , {<z}}) the 
greatest (resp. smallest) element in Vq. 

The permutation group Sq admits a natural left action on the partitions of [q]. Call On the orbit 
to which n belongs. These orbits are in natural one to one correspondence with the partitions of the 
integer q. For any permutation a € Sq, we denote by [a] the partition of [q] whose blocks are the orbits 
of a. 

The lattice Vq admits a subset of interest, of set of all non-crossing partitions. We denote it by 
NC{q). A partition is called non-crossing iff there exists no two different blocks X,Y of the partition 
and a,c G X, b,d G Y, with a < b < c < d. NC{q) is a sub-poset of Vq and it turns out to be a 
lattice as well. The operation A is the same for NC{q) and Vq, however the suprenium operation V 
is different in the non-crossing setting, as it can be seen by the example 11 = {{1}, {3}, {2, 4}} and 
n' = {{l,3},{2},{4}}. 

We consider the vector subspace Ap of functions C^'-^'-*-'^^'^^''' whose value is when Hi ^ 112 
does not hold. This vector space is endowed with a convolution operation * 

[/i*/2](ni,n2) ^ /i(ni,n3)/2(n3,n2). (7) 

n3 s.t. ni^n3<n2 

This turns Ap into an algebra with unit (5ni.n2- Fo^ each Hi, 112 G V such that Hi ^112, there 
exists a partition 11 such that the interval [Hi, 112] is isomorphic as a lattice to the lattice [Og,n]. If 
On — {pi, ■ ■ ■ ,Pk) then define the Mobius function Mob by 

Mob(ni,n2) = n ((-i)'-'c,_i)^' (8) 

i 

where 

1 f2i\ _ {2iy. 



i + l\ij {i + l)lil 

is the i-th Catalan number (later we shall see that Catalan numbers are a particular case of a much 
more general class of combinatorial family, called the Fuss-Catalan family). We also define the Zeta 
function C which is equal to one when Hi ^ 112 and zero otherwise: 

C(ni,n2) = i. (10) 

In this paper we will need the fact that the convolution product of the Mobius function Mob and the 
Zeta function C is equal to the Kronecker delta function (5(111,112) — 1 iff. Hi =112: 

Mob*C = '5. (11) 

Note that if we had considered the equation ([7]) in the general lattice of partitions and had left the 
definition of C unchanged, we should have modified the definition of Mob as follows. 



Mob(ni,n2) = n((-ir'(*-i)!r 
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We finish this section by mentioning the fohowing fact, of crucial importance in what follows: if 7 
is a permutation with one cycle (i.e. a cyclic permutation), the set Af := {a G Sq , \jcr^^ \ + |c| — \ j\ = 
(7 — 1} is in one to one correspondence with the set of non-crossing partitions NC{q). If 7 = (1, . . . , q), 
which we can always assume up to an overall conjugation of the set J\f, the correspondence is given by 
CT — )• [ct]. The proof of this result is elementary and belongs to combinatorial folklore. 

4.2 Weingarten calculus 

In this section, we recall a few facts about the Weingarten calculus, useful to evaluate averages with 
respect to the Haar measure on the unitary group. 

Definition 4.1. The unitary Weingarten function Wg{n,a) is a function of a dimension parameter 
n and of a permutation a in the symmetric group Sp . It is the inverse of the function a i— > under 
the convolution for the symmetric group (^a denotes the number of cycles of the permutation a). 

Notice that the function <j 1— !■ n^" is invertible for large n, as it behaves like nP5e as n — )■ 00, were p 
denotes the number of elements in the permutation group Sp. We refer to j41| for historical references 
and further details. We shall use the shorthand notation Wg((T) = Wg(n, a) when the dimension 
parameter n is obvious. 

The following theorem relates integrals with respect to the Haar measure on the unitary group 
U{n) and the Weingarten function Wg. (see for example [30]): 

Theorem 4.2. Let n he a positive integer and {ii, . . . ,ip), (i[, . . . ,ip), {ji, ■ ■ ■ ,jp), (ji,---,jp) be 
p-tuples of positive integers from {1, 2, . . . , n}. Then 

/ Ui,j, ■ ■ ■ U^^j^U.'^f^ ■ ■ ■ Ui^^y^ dU = 

JU{n) 

E5i,i' ■■■Sii' 5i,i' ...Sii> Wg(n,T(T~^). (12) 

If p ^ p' then 

[ U^in ■ ■ ■ ^W.C^- ■ ■ WJ^ dU = 0. (13) 

Ju{n) " " 

We are interested in the values of the Weingarten function in the limit n — 00. The following 
result encloses all the data we need for our computations about the asymptotics of the Wg function; 
see [30] for a proof. 

Theorem 4.3. For a permutation a € Sp, let Cycles((T) denote the set of cycles of a. Then 

Wg{n,a) = {-ir-*- H Wg(n,c)(l + 0(n-2)) (14) 

cGCyclcs((7) 

and 

Wg(n,(l,...,d)) = (-l)''-ic,_i n (15) 

where Cj = ./i^'l; .. is the i-th Catalan number. 

A shorthand for this theorem is the introduction of a Mobius function Mob on the symmetric group, 
invariant under conjugation and multiplicative over the cycles, satisfying for any permutation a £ Sp: 

Wg(n,a) =n-(P+l'^l)(Mob((T) + 0(n-2)). (16) 

where \a\ = p — fj^cr is the length of a, i.e. the minimal number of transpositions that multiply to a, 
and Mob was defined at equation ((S]). We refer to [H] for details about the function Mob. 

Next, we recall briefly the results of [JT] for the convenience of the reader and in order to make the 
paper self contained. 
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4.3 Axioms of unitary graphical calculus 



The purpose of the graphical calculus introduced in |31) is to yield an effective method to evaluate 
the expectation of random tensors with respect to the Haar measure on a unitary group. The tensors 
under consideration can be constructed from a few elementary tensors such as the Bell state, fixed 
kets and bras, and random unitary matrices. In graphical language, a tensor corresponds to a box, and 
an appropriate Hilbertian structure yields a correspondence between boxes and tensors. However, the 
calculus yielding expectations only relies on diagrammatic operations. 

Each box B is represented as a rectangle with decorations on its boundary. The decorations are 
either white or black, and belong to S{B) U S*{B). Figure 5(a) depicts an example of boxes and 
diagrams. 



M » 



(a) 



M 



X M 



(b) 



c A 


• 


• 

A(E)B 






c B 


■ □ 


I 



(d) 



Figure 5: Basic diagrams and axioms: (a) diagram for a general tensor M; (b) trace of a (1, l)-tensor 
(matrix) AI; (c) Scalar product {y \ M \ x); (d) tensor product of two diagrams 



It is possible to construct new boxes out of old ones by formal algebraic operations such as sums or 
products. We call diagram a picture consisting in boxes and wires according to the following rule: a 
wire may link a white decoration in S{B) to its black counterpart in S*{B). A diagram can be turned 
into a box by choosing an orientation and a starting point. 

Regarding the Hilbertian structure, wires correspond to tensor contractions. There exists an invo- 
lution for boxes and diagrams. It is anti-linear and it turns a decoration in S{B) into its counterpart 
in S*{B). Our conventions are close to those of [46j|47]. They should be familiar to the reader ac- 
quainted with existing graphical calculus of various types (planar algebra theory, Feynman diagrams 
theory, traced category theory). Our notations are designed to fit well to the problem of computing 
expectations, as shown in the next section. In Figure 5(b) 5(c) and 5(d) we depict the trace of a 
matrix, multiplication of tensors and the tensor product operation. For details, we refer to [31j . 



4.4 Planar expansion 

The main application of our calculus is to compute expectation of diagrams where some boxes represent 
random matrices (e.g. Haar distributed or Gaussian). For this, we need a concept of removal of boxes 
U and U. A removal r is a way to pair decorations of the U and U boxes appearing in a diagram. It 
therefore consists in a pairing a of the white decorations of U boxes with the white decorations of U 
boxes, together with a pairing /3 between the black decorations of U boxes and the black decorations 
of U boxes. Assuming that V contains p boxes of type U and that the boxes U (resp. U) are labeled 
from 1 to p, then r = (a, /3) where a, /? are permutations of Sp. 

Given a removal r S Rem(I?), we construct a new diagram Vr associated to r, which has the 
important property that it no longer contains boxes of type U or U. One starts by erasing the boxes U 
and U but keeps the decorations attached to them. Assuming that one has labeled the erased boxes U 
and U with integers from {1, . . . one connects all the (inner parts of the) white decorations of the 
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z-th erased U box with the corresponding (inner parts of the) white decorations of the Q;(j)-th erased 
U box. In a similar manner, one uses the permutation /3 to connect black decorations. 
In [31], we proved the following result: 

Theorem 4.4. The following holds true: 

r=(a, /3)gRom[/ (13) 



Main result: computing the moments of ps and a flow prob- 
lem 



5.1 A formula for the moments of ps 

This section contains one of the main results of the paper. Theorem 15. II 

We present a general method for computing the moments of the random density matrix ps obtained 
by partial tracing a random graph state (jl]) over the subsystems labeled by T. Here Iltracc — {S,T} 
denotes a partition of the total number n = 2m of subspaces. For each block d of Ilvortox, we define 
the sets Si — S D Ci and Ti — T D Ci. These sets define a partition of each vertex Vi and contain the 
information on which subsystems of the vertex Vi are being traced out. Let Ei^i — {k Cz Ci \ 31 Cz 
Ci, I > k s.t. (fc, I) £ E} be the set of bonds with both ends in Ci. Additionally, we have to take into 



account the bounds between the different blocks. For j ^ i , let E. 
E}. Notice that \E, 



{keC, 



31 G Cj s.t. 



= E, 



and that dsi^j — dsj^i ■ Since the edges of the graph are not oriented 
Each unitary block Cj is partitioned in two ways: by the partial 
tracing operation, Ci = SiUTi, and by the type of bonds it contains: Ci = UjEi^j. 

For the case depicted in Figure [U the vertex partition Ilvcrtox is made of 3 blocks: Ci {1}, 



we shall put, for i < j, Eij = Ei^j 



C2 — {2, 3} and C3 = {4, 5, 6}. We represent it in Figure 6(a) For the block C3 for example, one has 
^3 = {6}, T3 = {4, 5}, Es^i = 0, E3^2 = {4} and E^^s = {5, 6}. 



,¥3 



H2 Ti-s He 
(a) 



1*5+6) 



• • • 



PS 



T3 D iX 

2 3 6 



_ j_ 

N3 



• 




• • 















2 3 



U3 

4 5 6 



(b) 



Figure 6: A marginal ps of the graph state represented in Figure (TJ where Hilbert spaces 7^1,7^4 and 
■H5 have been traced out (a). The same marginal is represented using the graphical notation of ^31 in 
(b). 



The marginal we are considering in Figure 6(a) is ps — Tr^^^^^^gj^^ |\E')(^'|. In the diagrammatic 



picture of [31 , the density matrix ps is represented in Figure [6 (b)[ 

In the general case, the constant in front of the diagram is meant to normalize the Bell states 
appearing in ^ and it has a value of -^d[„] A^". Using the graphical Weingarten formula (Theorem 
14. 4p . we obtain the following theorem, which contains the main result of this section. 
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Theorem 5.1. The moments of a graph state marginal ps are given by the exact formula 

ai,...,ak, i=l i=l 

i3i,...,0keSp ^^j-j 

Proof. For each independent unitary matrix Ui, (dT^iVl^'l)'^"' counts the contribution of the loops 

corresponding to partial traces, (dg-iVl'^'l)^''''^ corresponds to the moment product, (rf^j.^TV''^"!)'' 

corresponds to loops created by Bell states inside Vi and, for each couple i < j, (dEi^N^^^^^^'^^^' 
represents the contribution of the loops coming from the Bell states between vertices Vi and Vj . □ 

The Weingarten functions appearing in equation (jl7p are intractable at fixed dimension N. In the 
following section, we shall investigate the marginals ps in the asymptotic regime iV — > cx) (keeping 
di,d2, ■ ■ ■ ,dn fixed). 

In order to compute the dominating term in equation ()17p . we look at the exponent of N for each 
term in the sum (we use the first order asymptotic for the Weingarten function (fT6|) and the fact that 
for all permutations a G Sp, #cr + \a\ = p)'. 

-^+jzm{p-\i-'c.A) + m{p-\a,\)\+ J2 \E^s\iP-w'|3,\)+i2p\Eu\+i2\^^\i-p-K'l3^ 

i=l l^i<j^k i=l i=l 

(18) 

Using, for aU i, \S^\ + \T,\ = |C,|, J^j^^ + 2\Eu\ = \Ci\ and the fact that \Ci\ = n, we conclude 
that the general term in equation ([T7| has the following asymptotic behavior: 

N-^--^ (jl Mob(a- + o(l) j , (19) 

where 

fc fc 

Fo,^=El^'ll^^'"'l + l^'ll"'l+ E \E^,\\l3-'l3,\+Y,\a\\a-'p,\. (20) 

i=l l^i<j^fc 1=1 

One has to minimize the function F^^p over ai, . . . , a^, /3i, . . . , G Sp in order to find the dominating 
term in equation (jl7p . The following simplification can be made at this point: 

Lemma 5.2. The minimum of the function F^.p is the same as the minimum of the function Fp 
defined by: 

k 

Fp^Y.\^,\\^~'p^\ + \^^\\p^\^ E \E^Jm'PJ\■ (21) 

i=l l^Kj^k 

Proof. Since \Ci\ = \Si\ + |Ti|, one may use \Si\ times the triangular inequality |7~^Q;i| + |a~^/3i| > 
|7~^/3i| and then \Ti\ times he triangular inequality \ai\ + |a~^/3i| ^ to reduce the minimization 
problem of (|20p to that of (pij) . The two problems have the same solution, since i^^(/?i, . . . , /3fc) ^ 
Fa,p{oii, ttfe, /?!,..., /3fe) for all ai and, choosing ai — Pi for all i, one has Fp — Fa^p. □ 

This time, Fp is a function which depends only on k permutations /3i, . . . , /3fc g 5p, hence minimizing 
Ffj should be easier than minimizing Fa^p. In the next section, we will present a complete solution to 
the minimization problems (EU)) and (pij) using a graph-theoretical method. 
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5.2 The minimization problem on permutations as a maximal flow problem 



We will now exhibit a connection between the minimization problems (1201) and (j2ip and a maximum 
flow problem on a network. For an introduction to the maximum flow on networks, see [32j . Chapter 
27. We introduce a network (V, w) with vertex set V, edge set £ and edge capacities w. The network 
is associated to the minimization problem for Fp ()21|) in the following way. The vertex set is given by 
V — {id,7,/3i, . . . ,/3fc}, with two distinguished vertices: the source s = id and the sink t = j. The 
edges in £ are oriented and they are of three types: 

f = {(id,ft); |T,| >0}U{(A,7) ; >0}U{(ft,/3,),(/3„A) ; \E^J\>0}. (22) 
The capacities of the edges are given by: 

zi;(id,A) = |r,;|>0 

zi;(A,7) = >0 (23) 
w{P,,(3,)^w{l3j,l3,)^\E,j\>0. 



The network corresponding to Figures 6(a) and |6(b)] is represented in Figure [71 




Figure 7: Network corresponding to the graph state marginal of Figure 6(a) 



Let us now consider the maximum flow problem on the network {V,£,w), with source s = id and 
sink t — "f. A flow in a network is a function / : V x V — > M with the following three properties: 

1. Capacity constraint: for all vertices u,v eV, f{u,v) ^ w{u,v); 

2. Skew symmetry: for all u,v gV, f{u,v) — ~f{v,u); 

3. Flow conservation: for all vertices u different from the source and the sink w G V \ {s, t}, 

J2f{u,v)^0. (24) 
vev 



The value of a flow / is the quantity 

\f\ = Y,fis,v). (25) 

In the maximum flow problem, we want to determine a flow of maximum value in the network 
(y,£,w). The Ford-Fulkcrson algorithm states that solutions to the maximum flow problem are 
obtained as sums of augmenting paths. One starts with an empty flow / = 0. An augmenting path 
is obtained by pushing a quantity x of flow through a path id — ^ /J^^ Pi^ (3^ — >■ 7, such 

that, for all edges (u, v) of the path, x ^ w{u, v). One can then update the flow function by adding x 
units of flow to each of /(id, /^^J, /(/^i^ , /^^j), . . . , (/3i, , 7). The weights are also updated by subtracting 
X from the capacities of the edges above. The new network, with updated capacities, is called the 
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residual network. One iterates this procedure until no more augmenting paths exists. The maximal 
flow that can be sent from s to t is denoted hy X = |/|. The final residual network is denoted by 
(Vros, i^^rcs, ■fi'ros)- Notc that frcs coutains only edges for which w^cs > 0. 

Let us analyze now the connection between a solution of the maximal flow problem and Fp. Each 
time a flow x is sent through a path id —>■ —>■ /^ij —>■••• /3i, — >■ 7, we write the triangular 
inequality 

^[lAJ + l/3n'/3.J + --- + |/?-_\AJ + |/3-M] ^x\j\=x{p-l). (26) 

Summing up all this inequalities for the augmenting paths that add up to the maximal flow X, we 
have that 

+ J2 ^(/3.,/3,)l/5~V,|. 

Moreover, we claim that equality can be achieved in Fp ^ X(p — 1). The residual network has 
the property that it contains no augmenting paths. This means that the source id and the sink 7 
are in different connected components of the graph (Vrcsi^rcs)- The equality case in i^^ X{p — 1) 
is obtained by setting — id for all permutations /3i in the connected component of the source id. 
Similarly, put /3i = "f for all vertices f3i in the connected component of the sink 7. Finally, impose that 
A = Pj for all permutations /?i , 13 j lying in a connected component different from the ones containing 
id of 7. In this manner, all extra terms in equation (|27p are zero, and one has indeed Fp = X(p — 1). 

In conclusion, we have showed that Fp ^ X{p — 1) always holds, where X is the maximum flow 
on the network (V,f,w). Of course, X does not depend to the solution (i.e. the choice of the 
augmenting paths leading to the maximum flow) to the maximum flow problem. Moreover, we have 
showed that for each solution equality can be attained by choosing permutations (3i in a way 
depending on the residual network (which depends on the solution J-). Let us now describe exactly, 
for a solution J- to the maximum flow problem, the set of fc-tuples {f3i, . . . ,/3/j) S Sp which saturate 
the inequality Fp ^ X{p— 1). 

Firstly, each time an augmenting path id — >■ — >■ /3i, 7 in is chosen, we 

use the triangular inequality (|26|) . This implies that all permutations /3i^ , /J^j , . . . , /3i, are geodesic 
permutations (i.e. elements oiSNcip) ^^id that they satisfy 

[AJ«;[AJ<---<[AJ, (28) 

where we note by [a] the cycle partition of a permutation. Actually, since we are dealing with geodesic 
permutations, it follows (see e.g. the last paragraph of section |4T]) . that these partitions are non- 
crossing. Hence, each augmenting path imposes two conditions on the permutations it contains: they 
should all be geodesic permutations and their associated non-crossing partitions should satisfy a linear 
ordering inequality (|28p . Moreover, the condition that the extra terms in equation (l?fl) should be zero 
imposes additional constraints on the permutations (except, of course, in the case where £ios is empty, 
and all connected components of the residual network are singletons). These conditions impose that 
all permutation lying in the same connected component of (Vios, i^^rcs) should be equal (so that the 
terms are all null). Hence, the permutations (3i achieving the minimum for Fp, for a solution 

J- of the maximum flow problem, should satisfy the following constraints (note that we are referring 
to the connected components of the flnal residual network (Vres, i^^rcs)): 



[Ai] ^ [A2] =^ • • • =^ [Ail: for all augmenting paths in the solution J^; 

Pi — id, for all vertices (3i in the connected component of id; 

Pi = 7, for all vertices Pi in the connected component of 7; 

Pi = Pj, for all vertices Pi,Pj in the same connected component. 



(29) 
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Two important comments should be made at this point. Using the properties of the augmenting 
flow, it is fairly easy to see that the above system of equations and inequalities is consistent. Moreover, 
all solutions of (|29p are geodesic permutations. We have seen that this already holds for permutations 
f3i which belong to at least one augmenting path. A permutation not belonging to any augmenting 
paths is connected in the residual network to either the source id or to the sink 7 (but not to both, 
since in that case id — > /3i — > 7 will be a valid augmenting path in the final residual network, which is 
impossible). In that case, /3i = id or /3i = 7, which are both geodesic permutations. 

Let Bjr be the set of solutions to (I^Hl) . It can be described by a poset Vj^, endowed with a partial 
order ^jf: 

= {/3i, . . . , e SNcipf I [A] [ft] whenever i^^j}. (30) 

Since there are more than one possible solutions to the maximum flow problem (all having the same 
value X), the general solution to the minimization problem (1211) is given by the union over possibilities: 

B= y B^. (31) 

J-" solution to the 
Max-Flow problem 

The solution of the minimization problem for the /3's being settled, let us move now to the initial 
minimization problem for over ai, . . . , a^, . . . , stated in equation (j20l) . We solve the 

problem completely, by characterizing all the permutations ai,...,afc which achieve the minimum. 
We shall consider three cases for a vertex i (which depend only on the initial network): 

(I) \Si\ > 0,\Ti\ > 0: the vertex (3i is connected to both the source and the sink in the network 
iV,£,w); 

(II) \Si\ = 0, \Ti\ > 0: the vertex (3i is connected only to the source in the network (V, £,«;); 

(III) |Ti| = 0, jS'il > 0: the vertex (3i is connected only to the sink in the network {V,£,w). 

Definition 5.3. A vertex Vi is said to be of type "T" if all the subsystems of Vi are traced out. A 
vertex Vj is called a type "S" vertex if none of the subsystems of Vj are traced out. 

In the first case, when going from the i^Q,_^-minimization problem to the ^^-minimization problem, 
we used both triangle inequalities |7~-^ai| + |a~"^/3i| > |7~"^/3i| and \ai\ + \a~^f3i\ > at least once. 
For a fixed vertex i and a fixed permutation the unique solution to the system 

|7-V,| + |a-V,| = |7-iA| 
\a.\ + \ai'f3,\ = m 

is ai = Pi. This follows from the following chain of inequalities: 

p-l + 2\a-^P,\ < |7-ia,| + \ai^ + \a,\ + ja-^AI = I7-IAI + = 1, (33) 

and hence |a~^/3i| = 0, for all vertices of type (I). 

In the case (II) above, since 15^1 = 0, only the triangle inequality \ai\ + |a~^/3i| = 1/3; | must be 
saturated. Hence, for a given geodesic permutation the set of permutations which saturate the 
triangle inequality is the geodesic set id — ai — >■ Pi. In other words, we have to consider all the 
geodesic permutations a; G Snc{p) such that [ai] ^ [/?;]. By considering all the possibilities for ai, 
one needs to compute the sum 

Mob(a-^A). (34) 

[ft] 

Using the fact that the function Mob is related to the Mobius function on the poset of non-crossing 
partitions, one can show (see [12], chapter 10 and equation (jlip ) that the above sum is non-zero only 
if the above sum is trivial, that is, only if /3i — id. Thus, permutations ai, f3i corresponding to a vertex 
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which is connected only to the source (i.e. a type "T" vertex) must satisfy ai = f3i = id. Similar 
ideas lead to the conclusion that all type "S" vertices (case (III) above) need to satisfy ai — f3i — 7. 
Moreover, we have shown, in all three cases (I)-(III) above, that the permutations ai, . . . , which 
achieve the minimum satisfy = for all i = 1, . . . , fc. This greatly simplifies equation (j20p . since 
the Mobius functions are trivial: Moh{a~^ /3i) = Mob(id) — 1. 

In conclusion, since some of the solutions in equation pip cancel out, we define, for each solution 
J- of the maximum flow problem, the set 

Bjr = {Pi,...,l3keBjr\/3,^ id for type "T" vertices and P-j = 7 for type "S" vertices}. (35) 

We introduce the modified (smaller) set of general solutions to the minimization problem: 

B^ U B^. (36) 

J-" solution to the 
Max-Flow problem 

We sum up the preceding discussion in the foUowing theorem, the main result of this section. 
Theorem 5.4. 

ETr(p-,) = iV-^(--i) Y: il(ds.f''~''''U{dT.)*''- n idE.,)*'''^'''i[dc:- (37) 

{l3i,....,l3k)eB i=l i=l l^'i<Kfc i=l 

In the simplified situation when the parameters are trivial, di — IVi, the above theorem admits the 
following corollary, which provides a simple combinatorial formula for the moments of a graph state 
marginal. Both combinatorial quantities below (the maximum fiow X and the cardinality of the set of 
solutions ifB) can be computed from the network associated to the marginal. 

Theorem 5.5. Consider a graph state \'^) with the property that the relative dimensions of its sub- 
systems are unity (di = ■ ■ ■ = dn = I). The average moments of a marginal ps — Tr^ |5')(^'| are given 
by the simple combinatorial formula 

ETr(p^) = 7V-^(P-i)(#B + o(l)), (38) 

where X is the maximum flow associated to the marginal and B is determined by the set of augmenting 
paths corresponding to the maximum flow problem. 

Theorem 15.51 is very convenient in the sense that it turns the problem of computing moments 
into the problem of counting the number of solutions of a maximum flow problem. This is a rather 
unexpected mathematical connection between quantum physics and networking theory. 

Remark 5.6. When one interchanges the elements of the partition defining the marginal 5 -O- T, all 
the objects in the above discussion are replaced by their duals. In the network, the source and the sink 
are exchanged; the value of the maximum flow does not change, and all the inequalities for the /3 's are 
reversed. This means that geodesic permutations ai and j3i are replaced by their Kreweras complement. 
Hence, the asymptotic moments do not change. This was expected, since the spectra ofTrj \'^){^\ and 
Tr jc differ only by null eigenvalues. 

In this paper, we are making heavy use of the Fubin-Study measure. It is the most natural one, 
as it ensures norm preservation of the states. However it is not the simplest one, and it would be 
computationally simpler to deal with random Ginibre ensembles instead of unitary matrices. The 
following observation shows that the results are asymptotically the same: 

Remark 5.7. The graphical calculus that we recalled in section^J^ has a Gaussian counterpart J^5| /, 
and it follows by rather direct inspection that the Theorems \5.4\ and \5.5\ would yield the same asymp- 
totic results if one replaced in our model unitary matrices by (properly normalized) non-hermitian 
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standard complex Gaussian random matrices. Non-hermitian standard complex Gaussian matrices are 
less natural (and in particular they are not norm preserving). However the graphical calculus is much 
simpler, as the Weingarten calculus has just to be replaced by Wick theorem (i.e. Feynman diagrams). 
The reason why one obtains the same result, is that for the Gaussian case, one does not have a 's and 
j3's but just 13 's and one ends up directly with the minimization problem for F{l3). The fact that the 
leading coefficients are the same follows from the non-crossing Mobius inversion formula 

We are not able to give physical reasons why the Gaussian and the unitary model have the same 
asymptotic behavior in general and it seems non-obvious to us without a direct computation of mo- 
ments. 



6 Applications 

In the last part of the paper, we study applications of the theory developed in the preceding sections. 
We start by introducing the Fuss-Catalan probability distributions which appear as limit eigenvalue 
distributions in some special cases. We then study particular graphs (stars graphs, cycles) and the 
density matrix ensembles one obtains by partial tracing the corresponding random pure states. Finally, 
in Section [6.61 we study marginals of graph-states which lead to new probability distributions. 



6.1 Free Poisson and Fuss-Catalan distributions 

In Random Matrix Theory, the free Poisson (or Marchenko-Pastur) distribution ttc^' describes asymp- 
totically the spectral density of the normalized random Wishart matrices W — GO* E A4n{'C), where 
G is a random rectangular matrix from the Ginibre ensemble of size N x cN. The Marchenko-Pastur 
distribution is also the limit spectral density of rescaled density matrices from the induced ensem- 
ble [26l[43] . If A denotes an eigenvalue of a normalized random Wishart density matrix W/ Tr W of 
size N the probability density of the variable x — NX is asymptotically given by the free Poisson 
probability measure with parameter c > 0, 



7r(i) = max(l - c, 0)So + \\+c-2^c,\^c^2^ dx , (39) 

where 1a denotes the indicator function of the set A. 

In free probability theory, this distribution is also called the free distribution, and it has a 
semigroup structure with respect to the additive free convolution of Voiculescu: ttc^^ ffl tt^^^ = ""c+d 

(see, e.g. [H]). Moreover, ttc^-* can be characterized by the fact that all its free cumulants are equal to 
c. 

The mean purity for the Marchenko-Pastur distribution reads J x"^ cZttc^^ (x) — (? -\- c. One can 
compute the mean entropy of this probability distribution (note that in this work, logarithms are 
considered in base 2, log = loga) 



= \ ~x\o^xd^\'\x) = l I. - -r^^.', (40) 

J l^T ifO<c<l. 

In the case c = 1 it is equal to —1/2, so the mean entropy of a random density matrix p of size N 
generated out of a square Ginibre matrix G G 7MAr(C) behaves asymptotically {N — >■ oo) as log A^ — 1/2. 

Next, we recall a few facts about a series of probability distributions that was discovered recently 
in relation to quantum group theory in [34] . and which generalize the free Poisson distribution p9p . 
These distributions depend on a real parameter s > and are most easily characterized by their 
moments. Four our purposes it is convenient to extend the meaning of the binomial notation: for an 
arbitrary real number a and a natural integer k we set 

^q;\ Q!(q; — 1) . . . (a — A: -I- 1) 

' — -. (41) 



k k\ 
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Theorem 6.1. For any real number s > 0, there exists a probability measure tt'^-*, called the Fuss- 
Catalan distribution of order s, whose moments are the generalized Fuss-Catalan numbers (see, e.g. 

FCl^^^^('P + P]. (42) 

The measure tt''*^ has no atoms, it is supported on [0, A'] where K = (s + l)'^^^ / s'^ , its density is 
analytic on {0,K), and bounded at x — K, with asymptotic behavior ^ l/(7ra;^/^'*+^-') at x — 0. 

This distribution arises in Random Matrix Theory when one studies the product of s independent 
random square Ginibre matrices, G = 11^=1 ^i- this case, s e N and squared singular values of G 
(i.e. eigenvalues GG*) have asymptotic distribution tt*^"). In terms of free probability theory, it is the 
free multiplicative convolution product of s copies of a Marchenko-Pastur distribution [34] : 



. (43) 



One can further generalize Fuss-Catalan distributions in the following way. Consider s independent 
rectangular Ginibre matrices Gi , . . . , Gg of dimension ratios ci , . . . , Cg and such that the product 
G — n^=i ^3 is well defined. The asymptotic level distribution for the (rescaled) eigenvalues of GG* 

is then ir^^ The new, generalized measure depends on s parameters c = {ci, . . . Cg} which describe the 
dimension ratios for each matrix. In free probability this corresponds to taking s Marchenko-Pastur 
distributions with different parameters c: 

7r(^)=i^W. (44) 

In the case of an integer s, the Fuss-Catalan numbers count the number of s-chains in the lattice 
of non-crossing partitions NG{p) [351: 

FC^P^ = |{6p < ^ ^2 ^ • • • < a, ip e NC{p)}\. (45) 

In [34], the authors interpreted these numbers as the moments oi free Bessel laws. They showed that 
these numbers count non-crossing partitions of [sp] into blocks of size multiple of s (for a bijective 
proof of this statement, see [33]): 

PCp ^ = |{^ e NG{sp) I V6 e T, |6| is a multiple of s}\. (46) 

Using the moment formulas, one can compute the entropy of the Fuss-Catalan distributions. 

Proposition 6.2. The entropy H of the Fuss-Catalan measure tt*-*-* is given by 

s+l 



F(7r('')) := j ~x\ogxdTT^'\x)=~^-. (47) 



Proof. Using the analyticity of the density, we can write the Shannon entropy as a limit of Renyi 
entropies, by replacing factorials in the definition of FGp by Gamma functions: 

/ log X d7r(-) (x) ^ hm log f^'l +^ + (48) 

J piii-p r(sp + 2)r(p + 1) 

Next, we write the above limit as a logarithmic derivative at p = 1 of the function 

T(sp + p+l) 

T{sp + 2)V{p+iy ^^^^ 
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Using the expression of the derivative of Gamma function at integer points 



r'(n + i) = n! l-^o + J2- 



(50) 



where 70 ~ 0.57721 is the Euler constant, we can conclude. 



□ 



The second member of the Fuss-Catalan family, tt*^^) corresponds to a product of s = 2 independent 
Ginibre matrices. Note that a related analysis of a product of two real random correlation matrices 
was performed by Bouchaud et al. [37j , who used the term Marchenko-Pastur square distribution, and 
independently by Benaych-Georges [44], who studied multiplicative convolution of free Poisson laws. 
The formulas in equation (j42p for s = 2 allow to explicitly solve for the density of this measure: 



^(27 + 3V81- 12a 



127r 



(27 + 3V81 - 12a;) 



1(0,27] (a;) dx. 



(51) 



A plot of the density of tt^^^^ is given in Figure [51 along with the density of the free Poisson 
distribution tt'"'^-' . 



density 




Figure 8: Plot of the densities for the probability measures tt*^^) and tt*^^). 



6.2 Marginals of graphs leading to Fuss-Catalan distributions 

We provide in this section examples of simple graphs and associated marginals whose limit distributions 
are members of the Fuss-Catalan family of probability laws introduced in the preceding section. The 
examples chosen here have the property of being the simplest ones which lead to Fuss-Catalan limits. 

The simplest case, s = 1, corresponds to the free Poisson (or Marchenko-Pastur) distribution. This 
distribution is the limit of the induced ensemble of random density matrices studied in [261143] . The 
simplest graph having the Marchenko-Pastur distribution as a limit, is the one-loop graph, with one 
vertex containing two subsystems (or Hilbert spaces) inside. The marginal we consider is the one where 
we trace out one system, keeping the other (since we assume that all the Hilbert spaces are isomorphic 
to C^, the indices of the systems being traced out do not matter). The simple graph, the marginal 
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^1 @ @ <7) 



(a) (b) (c) 

Figure 9: A vertex with one loop (a) and a marginal (b) having as a limit eigenvalue distribution the 
Marchenko-Pastur law tt'^^\ In the network (c), both edges have capacity one. 



and the associated network are represented in Figure [SI The maximum flow problem associated to the 
network is trivial, one unit of flow being sent through the path id — >■ /3i — >■ 7. The residual network is 
empty and it follows from Theorem 15 .41 that 



ETrp^^7Vi-f.|{[/3i]eiVC(p)|6p^[/3i]s;ip}| Vp^l. (52) 

The number of all non-crossing partitions of [p\ is known to be the p-th Catalan number FC^^ , which 
is also the p-th moment of the Marchenko-Pastur distribution it^^\ By construction both subsystems 
are of the same dimension so if one performs partial trace over one of them the resulting mixed 
state is distributed according to the Hilbert-Schmidt measure [42ll43] and the ratio parameterizing the 
Marchenko-Pastur distribution reads c = 1. 

We now move on to the flrst non-trivial case, s — 2. The graph we consider (see Figure [TU| has two 
vertices connected by an edge and each having a loop attached. The marginal of interest is the one 
obtained by partial tracing two copies of C''^ in the first vertex Vi and one copy of in V2 ■ In the 
network associated to this marginal, a maximum flow of 3 can be sent from the source id to the sink 7: 
one unit through each path id — >■ — >■ 7, i = 1, 2 and one unit through the path id — /3i — > /32 — 7- 
In this way, the residual network is empty and the only constraint on the geodesic permutations /3i , 
is 

Op ^ sc m ^ ip, (53) 
i.e. [/3i] and form a 2-chain in NC{p). It follows the moments of the marginal ps are given by 

and thus the rescaled random matrix N^ps converges in distribution to the second Fuss-Catalan 
measure tt'^^) . 



Vi V2 

(a) 




Vi V2 

(b) 

Figure 10: A graph (a) and a marginal (b) having as a limit eigenvalue distribution the Fuss-Catalan 
law TT^^^. In the network (c), non-labeled edges have capacity one. 

The construction for the graph state and its marginal in Figure fTOl can be easily generalized to obtain 
graph state marginals having tt^**^ as a limit distribution. The main idea behind this construction is 



the structure of the network in Figure 11(c) The maximum flow in this network is s -I- 1 and it can 
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be obtain by sending one unit of flow tlirough the following paths: 

id ^ ^1 ^ 7 
id ^ ^2 ^ 7 

(55) 

id I3s J 

id ^- /3i ^- /32 > -^1 

The residual network is empty, and thus the combinatorial factor appearing in equation (|17p is 
given by the number of s-chains in NC{p). In conclusion, the moments of ps behave asymptotically 
as 

ETrpP - Ar('^+i)(i-P'i^C^"). (56) 

The empirical eigenvalue distribution of the rescaled marginal N^^^ps converges to the Fuss-Catalan 
distribution vr'^'*^ . 




Figure 11: An example of a graph state (a) with a marginal (b) having as a limit eigenvalue distribution 
the s-th Fuss-Catalan probability measure 'k'^^\ The associated network (c) has a maximal flow of s -1-1, 
obtained by sending a unit of flow through each jSi and a unit through the path id — > /3i — > • • • — > 
Ps — 7- The linear chain condition [/3i] ^ • • • ^ [jSs] follows. 



6.3 One-unitary marginals of general graphs 

The next class of examples we shall investigate using the techniques developed in Section [3] is a rather 
general one. In this section, we consider marginals of graph states with the property that the set S 
of "surviving" subsystems is contained in one block of the vertex partition Ilvortcx- Note however that 
we do not impose any conditions on the underlying graph. 

We consider a general graph F with k nodes and m edges. As in Section [2J to such a graph we 
associate a random pure state \^) on a n = 2m- fold tensor product. We are interested in a marginal 

PS = TrT |*)(*| - Tr^,,,«, \^){^, (57) 

with the property that S* = [n] \ T C is contained in a single block h of the partition Ilvortcx defining 
the vertices of F. Without loss of generality, we can relabel the Hilbert spaces Hi, . . . ,Hti is such a 
way that b = {1, 2, . . . , n'}, where n' = \b\ is the number of subsystem of the "surviving" vertex. 

For a specific example, see Fig. [121 in which n' = 8. As we shall see, the other nodes of the graph 
F will not play any role in the statistical properties of ps, so they are not represented in Figure [T^ 

The single n'-vertex represents a random unitary matrix U, which acts on n' subspaces labeled by 
the set of indices in the block b. As before, T denotes the set of subspaces traced out, so the reduced 
matrix ps lives in the complementary subspace Hs- Let T' = [n'] \ S = T D [n'] the set of vertices in 
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b that are being traced out (T' C T). We are interested in computing the distribution of the random 
matrix ps- Since ps hves on Us C Hb, the probabihty distribution of the random density matrix ps 
is invariant with respect to any unitary conjugation on Us- 

Lemma 6.3. The distribution of ps is hi{'Hs)-i''^''Jariant. 




n_l n_l nj nj nt ni nj nj 

til rl2 rt3 iti 1X5 rte iLY rts 



Figure 12: Graph illustrating assumptions of Theorem 16.41 all but one vertices are completely traced 
out. The surviving vertex b has n' = 8 subsystems, out of which |G| = 4 lead to other vertices, while 
the remaining |F| = n' — \G\ nodes form closed loops. Out of these n' = 8 subsystems, \T'\ = 4 
represent subspaces which are traced out (crosses), while jS*] = 4 represent subspaces whose product 
supports the reduced state ps (full dots). The other vertices of the graph are not represented. 



Two subsets of [n'] = {1, . . . , n'} play an important role in what follows. The set S supporting the 
reduced density matrix and the set F of bonds which are contained in b: 



F= U {^,J}. 



(58) 



(ij)GBn[n'] 



In Fig. [121 S — {3,4,5,6} and F — {2,3,4,5}. In the graphical notation of [3T], the state ps 
corresponding to Fig. [T^] is represented in Figure [T31 The unitary blocks corresponding to vertices 
which were traced out were removed from the diagram, using the unitary axiom UU* = I. Using the 
box-manipulation rules of |31| . subsystems in b connected to other vertices in the initial graph are 
linked now by "identity wires" (subsystems 1,6,7,8 in Figure ITSl) . Notice that although the Hilbert 
spaces Hi with i S [n'] may have different dimensions, we used round-shaped labels for all of them, for 
obvious practical reasons. 

Let the set of bonds complementary to set F be denoted by G. It consists of \G\ nodes which 
represent bonds connecting nodes of c with nodes belonging to other vertices of F. In the case shown in 
Fig. [12] we have G = {1, 6, 7, 8}. Since the sets T' and S, and respectively F and G are complementary 
with respect to [n'], we have \T'\ + \S\ = \F\ + \G\ = n' 



• • • • 



PS 



T3 D D iX 

3 4 5 6 



u 

O O U 

1 2 3 4 5 6 7 8 




Figure 13: The diagram for the random state ps = TrT|5')(^'| obtained by reduction of the pure state 
shown in FigfT2l 

The normalizing constant 1/Z appearing in the diagram comes from the different normalizations 
of the input states: l/{diN) for every Bell state I'&J) associated to an "internal" bond {i,j) G F and 
l/{dkN) for every identity matrix corresponding to an index k G. Putting all these contributions 
together, we have 

Z=ll{d,Ny/'xll{d,N). (59) 
jgf jeG 
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The main result of this section is the following theorem, describing the asymptotic behavior of 
the marginal ps- We shall provide two proofs of this result: the first proof is direct, and it uses the 
graphical calculus in 31i. The second one makes use of Theorem 15.41 in Section [01 

Theorem 6.4. For a graph V , let |^) he the associated random pure state and consider 



PS = TrT|^')(*| 



(60) 



a marginal of this random state with the property that S is contained in a single block b of the partition 
^vertex defining the vertices of T . The vertex b contains n' nodes, out of which \G\ are connected to 
other vertices. The partial trace above is taken over \T\ subspaces of joint dimension N'^^'^dx- The 
reduced density matrix ps, supported by the complementary \S\ subspaces of dimension dsN'^^'^, is 
asymptotically (N oo) characterized by the following behavior 

(I) If \S\ < \T'\ + \G\, the state ps converges in moments to the maximally mixed state inHs- 

(II) If\S\ = |T'| + |G|, then the rescaled operator Ut'cIgN^^^ PS converges in moments to a free Poisson 
distribution tTc^"^ (see equation (j39p ) of parameter c = dT'dc/ds- 

(III) If \S\ > \T'\ + \G\, the reduced state ps has rank dT'dcN^^ /j ^/jg restriction of ps 

to its support, then ps converges in moments to the maximally mixed state in T-Lt' ® "Hg ■ The 
support of Ps is a dT'dcN^"^ dimensional Haar random subspace ofHs- 

Proof. Using the unitary invariance, it is easy to see that one can replace the tensor products of 
maximally entangled spaces by any rank-one projector. Using this trick and collecting all the vector 
spaces corresponding to the sets S, T' , F, G, one obtains the simplified diagram for ps depicted in 
Figure [Ml 



c Ps • 



z 



U 



U* 



Figure 14: Simplified diagram for reduced state ps corresponding to Fig. [T31 Round-shaped labels 
correspond to the subspace Hs of dimension dsN^^^, square-shaped labels to the traced out subspace 
Tix of dimension dxN^'^^, triangles to dpN^-^^ and diamonds to dcA^I'^l. The normalization constant 
reads Z = dcN^^^. 

Following Theorem 15.11 one can compute the moments of the random matrix ps using the Wein- 
garten function Wg(fc, a) and the graphical calculus: 



ETr(p^) = Z-^' [dsN^ 



s\ 



#(7-1q) 



#/3 



Wg(d[„,]7V"',a-i/3). (61) 



Since we are interested in the asymptotical regime N ^ oo, we investigate the power of N in the 
sum above: 

jV-p|Gl7vlS|(p-|7"'"l)7Vl'^'l(p-l"l)ArlG'l(p-|/'l)7vn'(-p-|""'^*l)_ (62) 

Using the equality n' = \S\ + \T'\, we find that the exponent of N that one wants to maximize in 
order to obtain the dominating terms is 



Using the triangular identities 



i\S\h-'a\ + \T'\\a\ + \Gm+n'\a~'p\). 

\r'a\ + \a-'l3\ ^ |7-i/3| 
|a|-f |a-i/3| > l/3| 



(63) 

(64) 
(65) 
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\S\ times and |r'| times respectively, we look at the (a priori) weaker minimization problem over /3 d Sp 

minimize /(/3) = |^||7-i/3| + (|r'| + |G|)|/3|. (66) 

Before solving this minimization problem, let us consider the extremal cases where 5 = or 
T' = 0. If S' = 0, then every subsystem of \'^) is traced out, and the reduced state is degenerate: 
P5 = 1 e A^i(C). If T' = then, for all N , ps is distributed as a Haar random projector of rank 
cLqN^'^^ in M.^gjyi\s\{C). From now on, we assume that \S\^ \T'\ ^ 1. We go back to the minimization 
problem (|66p . which is easily solved: 

min /(/3) = (p - 1) • min{|5|, \T'\ + |G|}. (67) 

PtOp 

The set of permutations /3 which reach the above minimum is as fohows: 

[{id} if|5|<|r| + |G|; 

argmin /(/?)= <^ {/? I id 7 geodesic} if |5| = |r'| + |G|; (68) 

[{7} if|5|>|r| + |G|. 



Having solved the minimization problem (j66|) . let us now return to the original maximization 
problem in a and /3, equation (|63p . Since (except in the trivial case where S or T' is the empty set) 
both triangular inequalities (|64p and (|65p have been used at least once, and since the only permutation 
a which saturates both inequalities at fixed /3 is a = /3 (see equations ([32| - ([33ll ). the dominating 
terms in the Weingarten sum are those for which a — j3 and /3 is optimal as in equation (|68p . We 
discuss now each case in equation (|68p separately: 



1. Case (I): \S\ < \T'\ + |G|. The dominating term is given by a = /3 = id. We get an expression 
for the expectation value of the power of the traces of the mixed state under consideration. 

ETr(p^) ~ 7Vl^l(i-P)d-Pd54'rf?;rfr«^] = (rfsA^l^l)'"' . (69) 
Thus, Ps behaves as a maximally mixed state on the dsiVl"^l-dimcnsional subspace Us- 

2. Case (II): \S\ — \T'\ + |G|. Here, dominating terms are indexed by geodesic permutations 
id a — f3 ^ "f. We get 



ETr(p^,)^iVl^l(i-^)d-^ Y: df''-'''Upd*-d-^^^ 

id— >-Q— >-7 



id— >-a— ^7 



(70) 



We conclude that the (renormalizcd) matrix dT'dcN^^^ ps converges in moments to a free Poisson 
distribution ttc^^ of parameter c = dT'dc/ds (see equation ([39|)). 



3. Case (HI): 151 > \T'\ + |G|. Here, the dominating term is the one with a — f3 = j. The 
asymptotic moments of ps are given by 

ETr(pP) ^ (dT'rfGiVl^'l+l«l)'"^ (71) 

We conclude that ps behaves like a maximally mixed state of dimension dx'dcN^^ Notice 
that this is the maximum rank of p^ in this case, since we are partial tracing over a space Ht' 
of dimension dx'N^'^ ' a state of rank dcN^'^V Hence, ps^ restricted to its support, converges 
in moments to a maximally mixed state. The distribution of the support of ps is Haar, by the 
unitary invariance Lemma [ 
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□ 



We shall now present a second proof of this result, using the flow network method of Theorems 
15.11 and 15.51 We invite the reader to judge the efficiency of the network method in finding the set of 
optimal permutations. 



Proof. The network associated to a one-unitary marginal is presented in Figure 15(a) Since vertices 
/32, ■ . ■ , /3fe are connected only to the source, one should have /32 = • • • = /?fc = id in the end. Hence, 
we can just ignore these vertices: they can only be used to send a maximum flow of \G\ from id to 



The updated, simplified reduced network is represented in Figure 15(b) 




\r\ + \G\ 



(a) 



(b) 



1^1 



Figure 15: Network associated to a one-unitary marginal (a) and its equivalent form (b) obtained by 
ignoring the vertices /32 , • . • , /3fc , which are traced out entirely. 

The maximum flow in the above network is X = minHS*!, \T'\ + \G\} and the residual network 
depends on which of \S\ and \T'\ + \G\ is greater, see Figure [TBI 



(a) 



T'\ + \G\ - \S\ 
id) — -(Pi 



(6) @ (g) 

\S\-\T'\-\G 
(c) @ @ <7 

Figure 16: Residual networks for one-unitary marginals, in the following cases: (a) \S\ < \T'\ + \G\ ; 
(b) |5| = |r| + |G|;(c) |5|>|r| + |G| 

In the first case, where \S\ < \T'\ + \G\, the permutation /3i is connected to the source by an edge of 
positive capacity |r'| + |G| — \S\ > and thus /3i = id. In this case, we conclude that all permutations 
ai, Pi are equal to the identity id. Hence, the asymptotic moments of the one-unitary marginal ps are 
given by plugging into equation (fT7)) the values of the optimal permutations: 



E 
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In the above formula, we have S = Si, T' = Ti and G — LJ*^^2^ij- After simphfying the factors in 
the products with the factors in Z, we are left with 



ETr 



i-p 



(73) 



which allows to conclude and is consistent with equation (|69)) . The third case, where \S\ > \T'\ + \G\ 
can be proved in an analogous way. 

Let us consider now the "critical" case, IS*] = |T'| + |G|. Since the residual network is empty (it 
has no edges), there is no constraint on the permutation /3i, beside that it is required to be geodesic. 
The permutations ai , /3i which contribute asymptotically are those such that 



ai = Pi 

Op < ^ ip 

oi'i = /^i = id Vi = 2, 3, , 



(74) 



Plugging these values into equation (jl7p . we recover the final result in equation (j70p . and the proof 
is complete. □ 



6.4 Star graphs 

In this section we investigate random quantum states associated to star graphs. A m-star graph is a 



graph with k = m + 1 vertices Vi, V2, . . . , Kn+i and edges {(j, m + (see Figure 17(a) ). For 

obvious graphical reasons, the vertices T^, V2, . . . , shall be called "satellites" and the distinguished 
vertex Vm+i shall be called the "center" of the graph. The central vertex is a tensor product of m 
subspaces 



Wm+l — ®i™m+l'Hi, 

while the satellites have only one subsystem: Wi — Hi, ior i = 1 



(75) 



, m. Edges correspond to entangled 
states between Hi and Hm+i for i = 1, . . . , m. We shall consider the simplified situation where di = 1, 
and thus Hi ~ for a ll i = 1, . . . , 2m. The graph state l^'star) G (C^)®^° associated to an m = 5- 
star is depicted in figure |17(b)[ Notice that the Hilbert spaces corresponding to the vertices are given 
by Wm+i = (C^)®™ and = for i = 1, 2 . . . , m. We have 




rn+l 



(g)|$+ 



171 -\- j ' 



(76) 



The unitary matrices Ui, . . . , Um+i are independent Haar-distributed random unitary matrices Ui G 
U{N), i = 1, . . . , m and U„,+i e W(iV™). 

Vi = ni 





Figure 17: A 5-star graph state in the simplified and usual graphical notation 

Before looking at general marginals of the star graph state, some simplifications can be made at 
this point. First, notice that one can assume Ui = ■ ■ ■ = Um = I. This follows from the fact that 
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the unitary transformation Ui (E) ■ ■ ■ <E) Um can be "absorbed" into Um+i', for a Bell state 1$+) and an 
unitary transformation [/, 

(C/0l)|$+)($+|(t/®I)* = (I(8)I7)|$+)($+|(I(8)I7)*, (77) 

which is a consequence of the well-known fact that the commutant of the group G of local unitaries 
G = {U®U\U(i U{N)} is spanned by the identity and the maximally entangled state: G' = 
Span{I, !$+)($+ 1}. 

Second, since the system is invariant with respect to permutations of the satellites l,...,m, a 
general marginal is specified by a coupe {s,t) of integers from [m], where m — s is the number of 
satellites that have been traced out and m — t is the number of traced particles in the central block 
0. The resulting marginal is pt^-*) g (^(^N-^^is+t) . ^ copies of have "survived" in the satellites and 
t copies in the central vertex. Let us make a remark about the rank of the matrix p^^'*\ Obviously, 
we have rkp'-"'*-' ^ s + t. Moreover, since we are partial tracing a pure state over (m — s) + (m — t) ~ 
2m — {s + t) copies of C^, we also have rkp^'*-*^ ^ 2m — {s + t). In conclusion, we obtain the bound 
rkp("^*) ^ min(s + t,2m- {s + t)). 

We shall look in detail at three examples, for m = 2-star graphs: {s = 0,t — 2), {s = l,t — 0) and 
{s = t = 1). The graphs corresponding to these marginals are represented in Figure fTSl We use the 
graphical formalism of ^Slj to represent the (random) density matrices p^^'^'' and p*-^'^^ in Figure 

m 




(a) (b) (c) 



Figure 18: Graphs corresponding to the marginals p^^'^^ and p'^^'^^ respectively. 



In Figure 19(a)[ we use the rules for manipulating boxes in the graphical formalism to show that 



p^"'^-* = Iat^ . Actually, we start by collapsing the Vi and the V2 boxes (and their adjoints) and then 
we use the same rule to collapse V3 and V3* . These considerations can be easily generalized to obtain 
the following lemma. 

Lemma 6.5. For all s,t G [m], the marginals pC^'*) and p^'^'^'> are given by the maximally mixed states 
in ■^In': G A^Art(C) and -^In^ G A^Ars(C) respectively. 

Let us analyze the marginal p^^'^\ The network corresponding to the minimization problem asso- 



ciated to this marginal is represented in Figure 20(a) 



Since the vertex (32 is connected only to the source, and the other two vertices /Sa and (3i are only 
connected to the sink, the solution for the minimization problem is given by 

fas = ^3 = ai = ^1 = 7, ^^g^ 
I a2 = = id . 

The maximum flow in the network is 1, and thus the asymptotic moments of the marginal are given 
by Theorem 15.51 

lETr(p(i'2))P^_l. + o(7V-P). (79) 

In other words, the rescaled random density matrix Np^^'^") converges in distribution to the Dirac mass 
Si. In physical terms, the random matrix p^^'"^^ (which has rank at most N) behaves asymptotically 
as a (normalized to unit trace) projector on a subspace of dimension N of . 

We finish our case studies by the marginal p*^^'^^ of a 2-star graph. The flow network associated to 
' ' is depicted in Figure [20(b)| The maximum flow in this network is 2: a unit of flow is sent through 
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Figure 19: Graphical notation for p^^'^^ and p^^'^^. 




Figure 20: Networks associated to the minimization problems for p(^'^)-(a) and 
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the augmenting paths id — )■ /Sa — >■ 7 and id — )■ /32 — >■ /^a — >■ /3i — > 7. The residual network is empty and 
thus, the geodesic permutations /3i,2,3 which achieve the minimum should satisfy /32 ^ ^ Since 
/3i is connected only to 7 and /32 is connected only to the source id, the sum over the Mobius functions 
constrain /3i and /32: 

ai = Pi = 7, 

a2 = /32 = id, (80) 
Op < [as = k] ^ ip. 

It follows that the asymptotic moments of the marginal p^^'^-* are given by 

>Tr(p(M))P = J_ J2 1 + (81) 

We conclude that the renormalized density matrix A^^p^^'^^ converges in distribution to a free Poisson 
distribution (of parameter c— 1) vr'"'^'. 

The above discussion can be easily generalized to obtain the following general theorem. 

Theorem 6.6. The marginal p^*^'*' g AAj^s+t (C) obtained by partial tracing m — s satellites and m — t 
central systems of a random star graph state {"^star) has the following asymptotic behavior: 

1. Ifs^O (resp. t^O), then, for all N, 

P^°^'^^^lm (resp. p(^'°^^^lMs). (82) 



2. If s + t < m, then iV'+'p*^"'*) converges in distribution to Si . 

3. If s + t > m, then p(*'*) has rank 2m — (s + t). If p^^'*^^ is the restriction of p^'^'^^ to its support, 
then iV^™~^*'+*'/5''''*-' converges in distribution to Si. 

4. If s + t ~ m and s,t ^ then iV^p'-"'*-' converges in distribution to a Free Poisson distribution 

Moreover, the average von Neumann entropy and the average purity of the marginals p(*'*) are given 
in the following table: 



Parameters 


von Neumann entropy E/f(p^*'*^) 


Purity ETr((p("^*))2) 


s = (resp. t = 0) 


ilogiV (resp. slogiV) 


N-*' (resp. N-') 


s + t < m 


(s + t) log iV + o(iV- 


AT-U+t) +o(iV-('^+*)) 


s + t > m 


(2m - s - t) \ogN + o(iV-(2'»-^-*)) 




s + t = m; s,t ^ 


mlogiV- 1/2 + 0(1) 


27V-" + o(iV-™) 



Tabic 1: von Neumann entropy and purity for marginals of star graph states. 



6.5 Cycle graphs 



In this section we look at random quantum pure states associated to cycle graphs. The asymptotic 
eigenvalue distributions of marginals of such graphs will be characterized in terms of the subset of 
traced systems. The set of possible limit measures one can obtain in this situation is fairly rich: for 
any classical multiplicative convolution of measure from the Fuss-Catalan family, one can construct a 
cycle graph having this measure as the limit eigenvalue distribution. 

An m-cycle is a graph having k = m vertices and m edges connecting the vertices in a cyclic way. 



The simplified and detailed representations of a m = 4-cycle are given in Figures 21(a) and 21(b) 
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Each vertex of the graph contains two subsystems, thus the whole graph has n = 2m subsystems. To 
keep notation simple, we shall assume that di — 1 for all i = l,...,2m. Hence, the vector spaces 
"Hi, ... , 'H2m are isomorphic to C^. Vertices are vector spaces 



= 'H2^-l ® n2^ =^ C 

The m edges of the graph correspond to maximally entangled states \^2j 2j+i) ^ 
(g) (with the obvious convention 2m + 1 = 1): 



(83) 



n2 



cycle/ 



(g)l* 



(84) 



where the random unitary operator Ui acts on Wi = 7^21-1 ^ 7^2 

Vi 




V2 




Figure 21: A 4-cycle graph state in the simplified and usual graphical notation 
We are interested in the eigenvalue statistics of the marginals 

PS = Ttt |*cyclo)(*cyclo|, (85) 

parameterized by partitions Iltracc — {S,T} of [n]. The number of systems being traced out in each 
vertex can be 0, 1 or 2, so we can classify the m vertices of a cycle in three classes: 

1. Vertices of type "S", where nothing is traced out; 

2. Vertices of type "R", where exactly one system is traced out; 

3. Vertices of type "T" , where both systems are traced out. 

Notice that we are implicitly using the fact that in a vertex with systems of equal dimension, the exact 
indices of the systems being traced out are not relevant; in this way, the definition of type "R" vertices 
is not ambiguous. 



In Figure 22(a) we consider a marginal of the m = 4-cycle, with S — {3, 4, 5, 8} and T — {1,2,6,7}. 



Vertex Vi is of type "T" , V2 is of type "S" and vertices V3 and V4 are of type "R" . We shall work out 
this example in detail, since it is the simples situation where the Fuss-Catalan distribution of order 2, 
tt'-^-' emerges. 

The network associated to the marginal ps is represented in Figure 22(b) The maximum fiow 
problem on the network has an unique solution which can be computed in the following way. First, 
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Figure 22: (a) Marginal of a 4-cycle graph with one "T" vertex (Vi), one "S" vertex (V2) and two 
"R" vertices {V3 and V4). (b) The network associated to the marginal; edges with no labels have unit 
capacity. 

send one unit of flow through the "R" vertices id ^ /^a ^ 7, id — > /34 — > 7. The only remaining 
way of sending flow from the source id to the sink 7 is through the arcs id — /3i — > /32 ^ 7 and 
id — >■ /3i — >■ /34 — >■ /Ja — > /32 -> 7- In this way, two additional units of flow can be sent from id to 7. 
After sending 4 units of flow (2 using "R" vertices and 2 using "TS" and "TRRS" arcs), the residual 
network is empty and the solution to the maximal flow problem is thus X — A. 

The set of geodesic permutations /3i , . . . , /34 which achieve the maximum flow is constrained only 
by the inequality [/3i] ^ [/34] ^ [/Ja] ^ [/32] (recall that for a permutation (3, [/3] denotes the partition 
induced by the cycle structure of /?). The fact that (3i (or V^i) was only connected to the source (i.e. 
it was completely traced out) imposes the additional constraint /3i = id. Similarly, since /32 is only 
connected to the sink (i.e. none of its subsystems was traced out), it must be that ^2 — 7- It follows 
that permutations , f3i which achieve the minimum in equation (|66p are those that verify 

Oil = Pi — id 

Using Theorem l5.51 we conclude that the moments of the above marginals are given by the following 
expression 

lim N^^P-^^ETr pP = |{cri,a2 € NC{p) \ cti CT2}|- (86) 

A''— ^00 

The combinatorial factor is the above equation is the number of chains of length 2 in NC{p), the 
Fuss-Catalan number FC'^\ see also equation p5|) . We conclude that the random density matrix 
N^ps converges in distribution to the Fuss-Catalan probability measure of order 2, 7r^^\ 

One can easily generalize the previous considerations to all m and arbitrary marginals to obtain 
the following the main result of this section, a complete characterization of the limiting eigenvalue 
statistics for marginals of random cycle graph states. 

Theorem 6.7. The asymptotic moments of the resettled random density marginal ps are products of 
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Fuss-Catalan numbers corresponding to "TR- ■ ■ RS" arcs: 

lim EiV('^«+l-^l)(P-i) Tiff, = TT FC(I"I) - n f^'"' ^ ^^^Y (87) 

where A is the set of all "TR- ■ ■ RS" arcs in the cycle graph, \a\ is the length of an arc a £ A and 
kn is the number of type "R" vertices. The empirical eigenvalue distribution of a rescaled version of 
ps converges to a classical multiplicative convolution product of Fuss- Catalan probability distributions 
Xag_47r(l'^l). Moreover, the average von Neumann entropy of the random density matrix ps is given by 

H+i 

EHips)^ikn + \A\)\ogN-J2 E -+«(!)• (88) 

Proof. For a general marginal ps of a random cycle graph state |^cycic)i let kji g ^ be the respective 
numbers of "R","S" and "T" vertices in the marginal; one has k^ + kg + k^ = k = m. Let also 
A be the set of all "TR---RS" arcs. An element a S .4 is a sequence of \a\ consecutive vertices 
^4o, V'oi, ■ • ■ , such that Va^ is an "T" type vertex, Va^^^^^ is an "S" type vertex and all inter- 

mediate elements , . . . , Va,^! are of type "R" (we consider only non-empty arcs, hence \a\ > 0). 
Arguing as we did in the previous example, the maximum flow in the general case is equal to kn -\- \ A\: 
one sends 1 unit of flow through each "R" vertex and 1 unit through each "TR- • • RS" arc a G A. The 
set of (geodesic) permutations {ai,/?^}™;^ which achieve the minimum in are characterized by 

at = l3t = id for all type "T" vertices Vt 
as = (is = 1 for all type "S" vertices Vs 

tti — Pi = id for all type "R" vertices Vi situated on a "TR- • • RT" arc 
aj ^ f3j ~ J for all type "R" vertices Vj situated on a "SR- • • RS" arc 
Op s=; Ki = l3aA ^--^ [aa|„| = I3a^j ^ Ip for all "TR- • • RS" arcs aeA 



The asymptotic moments of the marginal are obtained using our main Theorem 15.51 

^lim^iVC^^+l-^D^P-DETrp^ = J] l{l«|-chains in NC{p)}\ = [] FC^I'^D = [] ^ 

(89) 

The numbers in the right hand side of the equation above are the moments of the following product 
probability measure: 

TTs = x.e^^d'^l), (90) 

where x denotes the classical multiplicative convolution of measures. Recall that if two indepen- 
dent random variables X and Y have respective distributions p, and ly, then their product XY has 
distribution p x i/. □ 

Note that the above result depends only on the partition Xltraco = {S, T} and can be read directly in 
the graphical representation. The limiting distribution can be inferred from the arcs of type "TR- • • RS" 
of the cycle graph. 

6.6 Exotic graphs 

In the previous sections, we discussed graph states which give raise to the following asymptotic spectral 
distributions: 

1. Maximally mixed states (these correspond to 6% — tt*^")), 

2. Free Poisson (or Marchenko-Pastur) distribution tt^^-'. 
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3. Fuss-Catalan distribution tt''^^ s ^ 2. 



All the asymptotic measures we have encountered up to this point are members of the Fuss-Catalan 
family {7r'-*'}sgN. Moreover, in the preceding section on cycle graphs, we obtained limit eigenvalue dis- 
tributions which are classical multiplicative convolutions of Fuss-Catalan measures, see equation (j90p . 
In this section, we shall exhibit genuinely new limit asymptotic distributions arising from marginals of 
random graph states. 

The moments of distributions from the Fuss-Catalan family count the number of chains in the 
lattice of non-crossing partitions 33J , see also equation (H5|) : 

(91) 



J d^(^) (x) = |{ai, . . . , a, e NC{p) \ Op ^ ai ^ ■ ■ ■ ^ as ^ ip}\. 



In the lattice of non-crossing partitions NC{p), the Hasse diagram of such a chain is presented 
in Figure 23(a) for the case s = 3. Classical multiplicative convolutions of Fuss-Catalan measures 
correspond to multiple disjoint chains from Op and Ip, see Figure [23(b) [ for the Hasse diagram associated 
to TT^'^^ X Tr^^\ The simplest case which can not be described by the Fuss-Catalan statistics and 
independent products comes from the Hasse diagram in Figure [23(c)[ Next, we investigate graph state 
marginals corresponding to this diagram. 
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Figure 23: Examples of Hasse diagrams in NC{p): tt^'^^ tt^"^) x tt^^^ and an exotic distribution. 



An ad-hoc graph which yields such a Hasse diagram is depicted in Figure 24(a) The particular 
marginal of interest ps = Try | ^exotic) (5* exotic | is drawn in Figure [24(b)| with the associated network 



in Figure 24(c) 

The maximum flow problem for the network in Figure 24(c)| can be easily solved: one can send 5 



units of flow from the source id to the sink 7 using the following augmenting paths: 

• id /3i 7 

• id — 5- /32 — > 7 

• id ^ /33 ^ 7 

• id /3i ^- /32 -> 7 

• id ^ /3i ^3 7 

The residual network is empty and thus the set of permutations Ui^Pi achieving the minimum in 
(|66p is given by the conditions = Pi for i = 1,2,3 and by the inequalities of the Hasse diagram in 
Figure 23(c) with Ui = j3i. The next proposition follows from our main result. Theorem 15.51 
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g 9 g (^V^'^^^^V^'^^^^f^) 




Figure 24: Graph yielding an "exotic" limit distribution: (a) the simple form; (b) the marginal yielding 
the exotic limiting distribution; (c) the associated network. 

Proposition 6.8. The asymptotic moments of the marginal ps are given by 

lim l-ETT{N^psr = |{'^i,^2,CT3 e NC{p) \ ai ^ a^, ai ^ a^jl (92) 

N—^oo iV 

These moments are the moments of the following probability measure 

7rexot»c = ^^'^K(7rW x^«), (93) 
where M is the free multiplicative convolution and x is the classical multiplicative convolution. 

Notice that the probability measure in equation (j93p is genuinely new, and can not be obtain from 
the Fuss-Catalan family by taking classically independent products. In a similar manner, one can 
construct graph state marginals with eigenvalue distributions obtained from the free Poisson measures 
via arbitrary classical and free multiplicative convolutions. Analytical and combinatorial properties of 
such measures will be investigated in some future work. 



7 Concluding remarks 

In this work we have introduced ensembles of random graph states and established some of their 
basic properties. Any graph consisting of k vertices and m bonds represents an ensemble of pure 
quantum states, which describe a system containing 2m parts. Each bond represents a maximally 
entangled bi-partite pure state, while each vertex of the graph represents a random coupling between 
the corresponding subsystems. In the simplest case of the model all subsystems are assumed to be of 
the same dimension N, which can be treated as a free parameter of the model. Note that in contrast 
to generic, structureless random quantum states studied by Page [5], the random states analyzed in 
this work do posses certain topological structure determined by the graph selected. 

Dividing the entire multipartite system into two disjoint sets, one can analyze the typical cor- 
relations between these parts. Technically, one studies the entropy of entanglement between both 
subsystems, defined as the von Neumann entropy of the reduced density matrix. The key result of this 
work consists in developing efficient techniques which allow one to establish the average entropy for 
an ensemble of states associated with a given graph and its concrete partition. 
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In the limit of large N the spectra of random states obtained by this procedure can be classified. 
We have shown for which cases the spectrum can be described by the Marchenko-Pastur distribution 
7r'^^\ also called free Poisson distribution. This universal distribution describes statistics of the Wishart 
matrices W = GiGl, where Gi is a non-hermitian random Ginibre matrix. 

In certain graphs states the partial trace leads to density matrices with spectra described by the 
Fuss-Catalan distribution Tr^'*^ studied earlier in [34]. These distributions are characteristic of random 
matrices with the structure GG* , where G = Oi=i ^ product of s independent Ginibre matrices. 
Moreover, we identified other examples of the graphs, for which the spectra are described by an even 
wider class of "exotic" distributions. 

The model presented here can be generalized in many directions. Instead of a maximally entangled 
pure state |<i>+) = J2i=i ^7^1*) ^ I*) belonging to Hd <E) Hd, each edge of the graph could represent 
a generic random entangled state |$) = X^iLi V^l^) ® K)' such that the non-negative numbers pi 



sum to unity. For instance, the random vector p may be generated according to the Hilbert-Schmidt 
measure [42], which corresponds to taking a random state of size (P. Such a generalization makes 
the model more realistic, but it should not change qualitatively the results obtained, since the mean 
entanglement entropy of the corresponding bi-partite system changes from logd to logd— 1/2. 

Any graph analyzed in this work may also represent a physical system in a different way. In this 
'dual setup' any bond represents a random unitary matrix which couples two corresponding subsystems, 
while any vertex with b bonds represents a &-particle GHZ-like entangled state. These ensembles, 
characterized by bi-partite interaction, are closer related to various physical models of interacting 
spins. Analysis of spectral properties of random states defined in a the latter procedure is a subject of 
a following publication. Observe that for any cycle graph, in which each vertex has exactly two bonds, 
both definitions yield exactly the same ensembles of random pure states. 

We conclude the paper with a rather general remark. Although the approach presented in this work 
is not directly related to any concrete Hamiltonian model of the physical interaction, it is capable to 
describe generic coupling between the subsystems. Not knowing any details about the kind of the 
interaction one assumes therefore that it can be mimicked by a random unitary matrix. Averaging 
over the Haar measure on the unitary group we obtain rigorous bounds for the average entropy of 
entanglement between any two specified fragments of the system. In the limit of large system size 
the bounds derived become exact, and are characteristic of generalized ensembles of random Wishart 
matrices. The results obtained for the average entropy are universal in the sense that they depend on 
the topology of the graph and its partition, but do not depend on the kind of the interaction between 
subsystems. Establishing a direct link between the results obtained in this work and the properties 
of typical matrix product states analyzed for concrete physical models fISJ remains as a subject of 
further investigations. 
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